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Abstract 

We present analytic force-free solutions modeling rotating stars and black holes immersed 
in the magnetic field of a thin disk that terminates at an inner radius. The solutions are exact 
in flat spacetime and approximate in Kerr spacetime. The compact object produces a conical 
jet whose properties carry information about its nature. For example, the jet from a star is 
surrounded by a current sheet, while that of a black hole is smooth. We compute an effective 
resistance in each case and compare to the canonical values used in circuit models of energy 
extraction. These solutions illustrate all of the basic features of the Blandford-Znajek process 
for energy extraction and jet formation in a clean setting. 
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1 Introduction 


The extraction of the spin energy of a rotating compact object by its plasma magnetosphere is 
a basic physical process that is likely at work in a number of energetic astrophysical phenomena, 
from pulsar winds and jets to active galactic nuclei (AGN), and possibly including gamma-ray 
bursts and tidal disruption events. For pulsars, this process is essentially the original unipolar 
induction of Faraday pQ, in which a conductor rotating in a magnetic held generates an electric 
held, whose voltage will drive a current if wires are present. The neutron star is the conductor and 
magnet, while the plasma functions as wires. The study of this process began in the original work 
of Goldreich and Julian [2] and continues today. 

That such a process can also occur for black holes |3] is a priori rather shocking. After all, black 
holes cannot source a magnetic held, cannot provide charge carriers, and cannot support outhow of 
local energy. Nevertheless, a number of physical processes conspire to allow black holes to behave as 
unipolar inductors: an accretion disk can provide the magnetic held, immersing a spinning black hole 
in such a magnetic held produces electric helds (with nonzero E- B) |1], these electric helds can pair 
produce plasma from the vacuum by accelerating stray charges to energies above the pair-creation 
threshold [5], and spinning black holes can support outhow of the energy relevant to physical work 
at inhnity Blandford and Znajek (BZ) connected these formerly disparate ideas in their 

ground-breaking paper [3] and demonstrated an astrophysically viable Penrose process—the analog 
of unipolar induction—with the potential to power AGN jets. 

BZ, along with much subsequent work, considered the force-free approximation for the plasma, in 
which the matter degrees of freedom do not contribute to the dynamics, leaving a closed, nonlinear 
evolution system for the electromagnetic held alone 0[ini[n]. While this is likely a good approxi¬ 
mation for much of the magnetosphere (and essentially all of it in the pulsar case), it also has the 
virtue of being the simplest set of plasma equations that can support the basic mechanism of unipo¬ 
lar induction by stars and black holes in plasma. In particular, while nonlinear, the equations are 
simple enough to admit analytic solution in a number of cases. Indeed, in the last few years there 
has been an explosion of progress on analytic solutions, due in part to new techniques imported from 
the relativity [HI [121 1131 El El El El E] and high-energy theory [191 EQI 121] communities. While 
numerical simulation is essential to gaining enough detail to eventually capture all of the relevant 
effects in real pulsars and AGN, analytic work is invaluable for elucidating basic mechanisms and 
efficiently exploring their range of operation. 

In this paper, we continue the effort to push analytic solutions as far as they can be taken. We 
begin in flat spacetime, which is a decent approximation for neutron stars and, following the method 
of BZ, a springboard for approximate black hole solutions. For stationary, axisymmetric, energy- 
extracting, force-free magnetospheres, there are three previously known families of exact solutions, 
which may be designated by the shapes of their poloidal held line^ vertical, radial, and parabolic0 

'^Toroidal refers to the (p direction, while poloidal refers to the remaining directions: p and z in cylindrical 
coordinates. Poloidal field lines are simply the projection of the magnetic field lines onto the poloidal plane. 

^The radial solution was found by Michel [32] and the parabolic by Blandford [33|. Vertical solutions appear 
in the analysis of relativistic jets specialized to full cylindrical symmetry, where the theory becomes an ordinary 
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Figure 1: Poloidal field lines of known families of exact force-free magnetospheres, with current 
sheets depicted by a thick line. Our new, hyperbolic solution includes the radial and vertical 
solutions as limits. For the plots we take 0 = 0 for simplicity, but the rotating solutions are similar. 


In this paper we will add a fourth, “hyperbolic” solution, whose basic geometry resembles that of 
an AGN, with a thin accretion disk sourcing the magnetic field. The four field line geometries are 
illustrated in Fig. 

The new family of solutions is parametrized by the disk’s inner termination radius b as well as a 
free function 0('^), corresponding to the angular velocity of each field line V’- (These concepts will 
be defined precisely below.) We can produce a complete model of unipolar induction, including jet 
formation, by choosing G('0) to equal the angular velocity of a rigidly rotating conducting star (of 
radius R) on held lines that intersect the star and to vanish elsewhere. This generates a conical jet 
of opening angle = arcsin(i?/6) surrounded by a thin sheet of current howing in the reverse 
direction. This model could apply to a neutron star surrounded by a magnetized accretion disk or, 
as a basic mechanism, to protostellar jets. 

Following the method of BZ, we can promote this exact hat space solution to an approximate 
solution for a slowly spinning black hole. Aspects of this procedure were previously carried out by 
Beskin, Istomin, and Pariev [221127], who were not in possession of the exact hat solution. Taking 
the disk to terminate at the innermost stable circular orbit (ISCO) of the black hole, the opening 
angle is ~ 20°, and unlike in the star case, the return current is smoothly distributed in the jet, so 
that there is no current sheet. One could in principle utilize this qualitative diherence to determine 
the nature of the central object from observations of the jetj^ Figure [^illustrates these results. 

There has been much discussion over the years of the extent to which black holes behave like 
conducting bodies. In possession of nontrivial analytic solutions with the same magnetic geometry 
in each case, we are able to make a direct comparison. In particular, we can test the idea that 
the black hole has an effective resistance equal to that of free space (equal to 1 in our units, or 

differential equation that can be directly integrated (e.g., Refs. [24l [25] 1. 

^Note that the jet will only be conical as far away as the force-free approximation holds; beyond that limit, gas 
confinement mechanisms will likely hx the shape. 
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(a) star 


(b) black hole 


Figure 2: Jet formation via unipolar induction by rotating compact objects in force-free plasma 
with magnetic field provided by a thin disk. We provide exact solutions for stars in flat spacetime 
and perturbative solutions for slowly rotating black holes. The main qualitative difference is that 
the jet produced by a star is surrounded by a current sheet (shown in red), while the jet produced 
by a black hole is not. 

approximately 377 Ohrns)]^ In a limit where the two solutions can be precisely compared, we hnd 
that the effective total resistance in the star circuit is equal to 0.5, while the effective total resistance 
in the black hole circuit is approximately 0.77. Thus, the black hole does in a sense provide extra 
resistance, but it adds a much smaller amount than that of free space. 

These solutions also possess pedagogical value. While the main features of the BZ mechanism—an 
accretion disk sourcing a magnetic held that threads the black hole, producing a jet—are by now 
well-established, there has not previously existed a complete analytic treatment demonstrating all 
the features in a single setting. Our treatment also contextualizes the BZ mechanism within the 
more established physics of unipolar induction. We hope that the calculations presented here will 
help newcomers understand this beautiful process. 

The remainder of the paper is organized as follows. In Sec. we derive the new family of exact 
hat force-free magnetospheres. In Sec. we model unipolar induction by stars with disks, and in 
Sec. 1^ we hnd approximate solutions in Kerr and model unipolar induction by black holes with 
disks (i.e., the BZ mechanism). We conclude with a comparison of the star and black hole cases in 
Sec. Our metric signature for hat spacetime is (—,-b,-b,-b) and we work in Heaviside-Lorentz 
units with G = c = 1. 

"^This notion holds precisely within the membrane paradigm imiiaES] as a surface resistivity dictating the 
behavior of ficticious electric and magnetic fields (defined by using infinitely accelerated observers with suitable 
renormalization). The idea of assigning the whole black hole a true resistivity of 1 has later been applied in models 
for binaries |3IJ 132] ■ 
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2 Hyperbolic Magnetosphere 


We consider force-free field configurations, i.e., closed 2-forms F (V[aFhc] = 0) satisfying the force- 
free condition 


FabJ'’ = 0 , ( 2 . 1 ) 

where = VcF^^ is the charge-current density four-vector]^ This is the statement of vanishing 
Lorentz four-force density, pE + JxB = 0 and J-E = Q. We restrict our attention to stationary and 
axisymmetric solutions with nonzero poloidal magnetic held, which may always be written mm 
as 

27 rF = d'l/j A [d(j) — dt] -I— dz A dp, (2.2) 

P 

where (t, z, p, 0) are cylindrical coordinates and we use the orientation dt A d^; A dp A d^j^ Here, 
'0(p, z) is the magnetic hux function, dehned to be the upwarcQ magnetic hux through the loop of 
revolution (about the symmetry axis) at the point (p, z) in the poloidal plane. The level sets of 'ip 
are the poloidal projections of the magnetic held lines, or poloidal held lines for short; these lines 
are also equipotential surfaces. The quantity / is the total electric current howing upward through 
the (p, z) loop. That I is only a function of ip corresponds to the statement that the poloidal current 
hows along poloidal magnetic held lines. The toroidal magnetic held is also proportional to I by 
= I/{2np). Finally, VLpip) represents the angular velocity of each held line 

The energy howing upward through a loop labeled by pj is given by 




The diherential voltage drop across a held line is dV = Vtpip) dp). Thus we have dP ~ I dV at a 
diherential level, which is reminiscent of the circuit equation P = IV. However, here I is the total 
current enclosed, while dP and dV are the local energy hux and voltage drop rates. 


The form of Eq. (2.2) ensures that three of the four components of the force-free condition (2.1) are 
satished. The remaining component in hat space is the so-called stream equation for p) 


(1 - p^n^) v^p) --dpp)- p^ivpjfnn' + ir = o, (2.4) 

p 

^Reviews of force-free electrodynamics, including the machinery used here for stationary, axisymmetric solutions, 
may be found in Refs. [27] (primarily the 3-1-1 approach) and [11] (primarily the covariant approach). 

®The electric and magnetic fields are given by 2'kE = —HVp and 27ri3 = p~^{Wp x cj) + Ip)- Our conventions 
differ from Ref. [TT] (GJ) hy p = 2'!rpGj and from Ref. [27] (B) by ^ I = — vTtt/b, E = i^B/'s/dF, and 

B = where the VTr accounts for the change to Gaussian units. 

^In this paper, by “upward” we mean using a surface that intersects the symmetry axis in the northern hemisphere 
and has orientation induced by the right-hand rule on the bounding loop. 

®One physical consequence of this notion is that collisionless particles making up the plasma will be stuck to the 
rotating line. It is also convenient because field lines that intersect rotating conductors must corotate (see Sec. |^. 
See (e.g.) Ref. m for further discussion of the notion of field line angular velocity. 
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where a prime denotes a derivative with respect to V’- If dehnite expressions for and 
are chosen, Eq. (2.4) becomes a second-order, nonlineaij^ partial differential equation for -0. For 
analytic solutions, this is not a particularly promising approach, since nonlinear partial differential 
equations are generally difficult to solve. Furthermore, it is hard to see how to choose physically 
interesting /('0) and f2('0) before determining -0 itself. 


An alternative approach consists of introducing an auxiliary variable u{p, z) and assuming that 0, 
I, and are only functions of u, 


'ip = 'ip{u), I = I{u), Q = Q{u). (2.5) 

This makes the level sets of 0 the same as those of u, so that by picking u one chooses the shape 
of the poloidal held lines but not their density. There is no guarantee that a force-free solution 
corresponding to a particular shape will exist, but at least one can choose physically interesting 
shapes and hope for the best. This is the method used by Blandford [23] to hnd the parabolic 
solution and by Menon and Dermer [3S| to hnd their null solution in Kerr. The radial and vertical 
solutions can also be derived straightforwardly this way. 

Our choice of m(p, z) is derived from the confocal ellipsoidal coordinate system for hat space. The 
coordinate lines are ellipses and hyperbolae with foci at a distance b from the origin. The curves of 
interest to us are the hyperbolae that arise as level sets of 


u = 


+ - \/(P - + 2^2 


26 


( 2 . 6 ) 


which have a shape that could plausibly resemble the held lines in an AGN. Figure plots these 
hyperbolae along with the corresponding value of u G [0,1]. Plugging Eqs. (2.5) and ( 2.6[ ) into 
Eq. (2.4), we are able to satisfy the equation provided that 


0'(m) = 00 


u 


_^00 jl + h‘^c[l-h‘^Vt{uYu‘^ 


(2.7) 


where 0o and c are constants and VL{u) is an arbitrary function. The choice of sign corresponds to 
the direction of current and energy how. In order to avoid a line current on the symmetry axis, 
I{u) must vanish there, J(0) = 0. This requires c = —6“^, leading to 


0 = 00 


u du 




I = ±00 


Q{u)u‘^ 


a/1 — 62G(m)2m^’ 


( 2 . 8 ) 


where the constant in the integral is to be hxed by requiring that 0(0) = 0 to make 0 the hux 
function. For G = 0, we can perform the integral explicitly, resulting in 

0|t^=o = V'o (l - . (2.9) 

^Choices of I and ft that linearize the equation have various physical difficulties (e.g., Sec. 2.6.2 of Ref. [27]i. 
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Figure 3: Level sets of m(p, z). The coordinate u ranges from u = 0 on the axis 6* = 0 to u = 1 on 
the disk 6* = 7r/2, p > h. Each line u ^ 1 meets the equator z = 0 ai p = bu. 


This vacuum solution of Maxwell’s equations was found previously in Ref. m- For a nonzero 
constant we can represent the integral with elliptic functions (see Appendix]^, while for arbitrary 
Q{u) we can proceed no further. However, the electric and magnetic helds are directly expressible 
in terms of and therefore we have an explicit closed-form solution, 


F = 


(V’o/27r) 




u du 


u 




A [d0 — f2(M) dt] ± f2(n)— dz A dp 


P 


( 2 . 10 ) 


Note that u{p, z) depends on b via (2.6). Equation ( |2.10[ ), or equivalently Eq. (2.8), provides a 
new exact family of force-free magnetospheres, parametrized by disk termination radius b, held line 
angular velocity f2(M), and current/energy how direction ±n 


2.1 Properties 

In order for the prefactor [1 — b‘^fl{uYu‘^] in the solution to be hnite and real, each held line 
u must satisfy bu^ < Rl{u), where Rl{u) = is the cylindrical radius of the light surface. 

Since p = bu is the place where the held intersects the equator and n < 1, the condition is satished 
if all held lines rotate subluminally on the equator, and in particular if they are anchored there to 
a rotating conductor (see the discussion in Sec. |^. 

The apparent singularity at the disk u = 1 is just a discontinuity in physical coordinates (e.g., 
cylindrical), and corresponds to the presence of surface charge and current. In particular, the 
pullback of F to the disk is continuous, indicating the absence of magnetic monopoles, but the 

^°One way to compute the electric and magnetic fields associated with the solution is to begin with ( |2.8[ ) and use 
the formulas given in footnote 
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pullback of the dual -kF is discontinuous, revealing the presence of surface charge and current. 
Using the formulas (A23) and (A24) from Ref. [H], we compute the surface charge and current to 
be 


Jsmi = J5{z)H{p — 6) dz A dp A (dt — hli dc/)) , 


J = 




( 2 . 11 ) 


where 6{x) is the Dirac delta and H{x) denotes the Heaviside function (vanishing for x < 0 and 
unity for x > 0). This corresponds to a surface current density Jsurf and charge density psurf of 


"Aurf iT<^('2)hf(p 0; Psurf pf^l 




surf 


( 2 . 12 ) 


Here, Di = is the value of D as the field lines collapse the disk, which would normally be zero 

in physical models (e.g., those of Secs. [^and|^. In this case, there is no surface charge, and the 
surface current is given by 


"Aurf 


(V’oA) 

Py/p2 - 62 


6{z)H{p-b)4>. 


(2.13) 


The current density blows up at the disk edge p = b (from the right), but the blowup grows only 
like the square root of the distance, so the total current is finite. Orthonormal frame components of 
the magnetic held also blow up like the square root of distance. This divergence is milder than the 
inverse-distance singularity that would be produced by a line current. Thus, our model contains a 
concentration of current near the edge of the disk, but not as much concentration as a line current. 
Similar comments apply to the charge and current density when Di ^ 0. 


Using Eq. (2.3), the power (energy hux) howing upward between two held lines ui < U 2 is given by 

du 


+2ir 


ru2 


’Ui 


a/1 — [1 — 62D(m)^M^] 


(2.14) 


We see that if the lower sign is chosen, the solution represents energy howing in from inhnity in 
the southern hemisphere and out to inhnity in the northern, while the upper entails the opposite. 
For physical models of unipolar induction (Secs. [^and|^, we have the freedom to choose the sign 
separately in each hemisphere for held lines that intersect the object and will do so in order to 
ensure outgoing hux everywhere. 

In spherical coordinates]^ the asymptotic r —)■ cx) behavior of the function u dehned in Eq. (2.6) is 


b^ 

u = sinO + 0\ -t: 


Thus, from Eq. (2.8), the large-r behavior of the solution is 

D sin^ 9 


I = ±V^o 


deij = i/jo 


\/r 


— sin'^ 9 
sin 9 


+ o|^ 


+ 0 ^ 


yi - iWsin'fl V f 


(2.15) 

(2.16) 
(2.17) 


Whe transformation to spherical coordinates is defined by p = rsin0 and z = rcosd. 
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Since 'ip is asymptotically r-independent, the field lines are asymptotically radial. They are uniformly 
distributed (monopolar) when = 0 and become more bunched near the equator as this quantity 
increases. From Eqs. ( |2.16| ) and ( |2.17[ ), we see that 


I = ±0, sin 6 deip, 


r —)■ cx). 


(2.18) 


This indicates that the electric and magnetic fields are asymptotically equal in magnitude. 

The solution includes the radial and vertical solutions as special cases. The radial solution emerges 
as 6 —)■ 0 at fixed spherical coordinates r and 6. From Eq. (2.15), we have u = sin 6 ' + 0{b‘^), and 
hence from Eq. (2.10), we find 

y ^0 . nf COS0 

hm p = — sm y - -— 

b^o 27r y I cosy I 


d^ A [ d0 — f2(sin 0) df] ± f2(sin 0) dr A d0 


(2.19) 


This is Michel’s radial solution |22], already in “split” form on account of the factor of cos0/ |cos 6 '|. 
Normally, one derives the true monopole solution and then adds such a factor by hand. Here, we 
see the factor naturally emerge when the solution is regarded as a limiting case of the hyperbolic 
solution. Note that if one chooses ± oppositely in each hemisphere so as to have pure outgoing 
(or pure ingoing) flux, then the solution becomes proportional to sin 6 * d 6 * A [ d0 — f2( df ± dr)] with 
opposite signs in opposite hemispheres, the simple form discussed in Ref. m- 


The vertical solution emerges as 6 —?• cx) at fixed cylindrical coordinate. 
u = p/h + 0{b~^), and hence from Eq. (|2.10 ), 


From Eq. (2.6), we have 


lim (b‘^F) = dp A [ d0 — H( df ± dz)] 

b^oo ' 27r 


( 2 . 20 ) 


In this particular limit, H is a constant because fl{u ~ p/b) -A- H(0), but one may choose H to scale 
with b so that the more general solution, Eq. ( |2.20[ ) with H = H(p), is obtained. Time-dependent, 
nonaxisymmeric generalizations of vertical solutions are given in Ref. ini. 


3 Star with a Disk 

Unipolar induction (also homopolar induction) is the generation of an electric field by a conductor 
moving in a magnetic field. The effect may be understood directly from the perfect conductor 
assumption that the electric field must vanish in the rest frame of the conductor. If the conductor 
moves with velocity v (and boost factor 7 ) in a magnetic field B, the rest-frame electric field is 
given by 'j{E + v x B), so that an electric field of —v x B must be generated. If wires (or plasma) 
are connected to the conductor, then the resulting voltage will drive a current, and energy will be 
carried away. In the case of an ordinary circuit, the amount of current and power is fixed by the 
resistance of the load. In the case of a conductor in plasma, one must self-consistently solve the 
equations, from which an effective resistance can be determined. 

In the framework of stationary, axisymmetric force-free magnetospheres, the boundary condition of 
a rotating conductor is simple and intuitive: field lines meeting the conductor must rotate with it. 
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This may be seen by recalling that the electric field in a frame with four-velocity m" is given by 
Fapu^. For a conductor rotating with angular velocity u, the four-velocity (field) is u" = 7 ( 1 , 0, 0,a;) 
in cylindrical coordinates. Plugging into Eq. ( 2.2[ ), we see that Fabv}’ = 0 requires that ^ 

for any field line "0 that intersects the conducting surface]^ This leads to the intuitive picture of 
field lines being anchored to the conductor. 


In light of the preceding discussion, it is straightforward to use our solution to model unipolar 
induction by a star of radius R < b rotating at constant angular velocity hi* in the field of a 
magnetized disk. We use the force-free solution outside the star and demand that field lines meeting 
the star corotate with it, 


^( m ) = n^H{R-bu), 


(3.1) 


where H{x) again denotes the Heaviside function. Second, we must impose the boundary condition 
that no energy comes in from infinity, by selecting the lower sign in the northern hemisphere and 
the upper sign in the southern hemisphere, ± —)■ —sign( 6 * — vr/2) in Eq. (2.10) [and Eq. (2.8) . 
Thus, the self-consistent model of unipolar induction by a rotating star in plasma in a magnetic 
field provided by a thin disk is 


F = 


where 


(V’o/2vr) 


\/l — b‘^Q{uyu'^ 


u di/, \ 

A [ d0 — fl^H{R — bu) dt] — s [R — bu) — dp A d^; ] , (3.2) 

H p J 


TT 


s = signl^ - - ) . 


(3.3) 


The model is symmetric under reflection about the equatorial plane, and for simplicity, we shall 

< 7 r/ 2 . The flux 

0o(l - \/l - u^), R <^b. 


henceforth consider solely the northern hemisphere, 0 < 6 ^ < 7 r/ 2 . The flux function is given bj-^^ 

u dw 


^(1 - ^ 2 ) [1 - b^niu^iR - bu)] 

Thus the total magnetic flux on the star is 

^ udu 

= 00 


(3.4) 


i>o f R 


’0 [1 - b‘^^lu^H{R - bu)] 


i? < 6 . 


(3.5) 


li R ^b, these integrals are given by elliptic functions (Appendix [A| . 

general, only the tangential components of the electric field must be continuous across a surface layer, whereas 
here we have demanded that all components be continuous (and hence vanish). If one considers only the tangential 
components, then 17 = w is still required provided that the field line intersects the conductor nontangentially. 
Furthermore, this makes all components vanish, so that there is no induced surface charge in the conductor frame. 
For further discussion, see Sec. 8.1 of Ref. m- 

Since the star cannot rotate faster than the speed of light, we have i?f7* < 1 and hence foPrt < 1, making 
<C 1 when R b. 
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3.1 Star Jet Properties 


The model contains current and energy flux only within the last held line touching the star, u < 
R/h. Since the held lines are asymptotically radial, this may be called a conical jet. In spherical 
coordinates, we have u ~ sin^ for large r from Eq. (2.15), so the (half-)opening angle is just 

R R 

0*. = arcsin — — R <^b. 

h h ^ 


(3.6) 


From Eqs. (2.14) and (2.15), the outgoing power per unit angle 9 in the asymptotic jet is 

sin^ 9 


dP 

~d9 


^_ 

27r 1 — sin'^ 9 ’ 


r —)■ cx). 


(3.7) 


which is concentrated toward the edge of the cone. (This formula includes only the contribution 
from one hemisphere; the total power is a factor of 2 larger.) Recall that Rl{u) = and hence 

612* = 6/Rl is the ratio of the disk radius to the light cylinder radius. 


Using Eq. (3.7), we can compute the total power from a slowly rotating star (612* 1), 


P* = 




2tt 


sin'^ 9 = 


Stt 


R2 


62 




1-TW 1+X7X 


262 J 


(3.8) 


If in addition R b, then 


P* 

-'jet 


8tt 


2ti 


(3.9) 


3.2 Circuit Model 


The net current flowing in the bulk of the jet is given by I{u = R/b), 


7buik — —'00 


12*(R/6)' 


Vl - (R12*)2(R/6) 


2 ’ 


(3.10) 


while an equal and opposite remainder flows along the current sheet at u = R/b. [This “return 
flow” is required by current conservation and may also be computed directly from the jump in 
caused by the discontinuous choice of 12 (m).] The voltage drop from the pole to the equator (or 
equivalently, to the end of the jet) is just 12*0and hence can be expressed in terms of elliptic 
functions (see Appendix]^. To leading order in 612*, it is simply 


V 




00 

27r 


yrqiT^l + . 


(3.11) 


This is the approximation that the light cylinder radius 1/12* is well outside the disk edge 6. Since 
R < b, this approximation also entails i?12* -C 1 small as well, i.e., that the surface velocity of 
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the star is much less than the speed of light. Taking this into account, we have for the effective 
resistance Res of the plasma in the circuit 


y - vi - (R/by 

/buikl ^ {R/by 


foff* -C 1. 


(3.12) 


Here, Res varies from 1/2 to 1 as R/b varies from 0 to 1. The value Res = 1 is the “impedance 
of free space” (dvr in cgs-Gaussian units), often used as the effective resistance of the plasma for 
rough estimates (e.g.. Ref. [S2])- Thus, for our setup in the slow rotation approximation, the true 
effective resistivity of the plasma differs from that of “free space” by at most a factor of 2. For 
rapidly rotating stars, one may use the expressions in Appendix The effective resistance never 
exceeds 1 and changes order of magnitude only when both fl^R and R/b are near 1, where it drops 
to zero. This is the limit where the star nearly touches the disk and has surface velocity nearly 
equal to the speed of light. 

If we instead assume R/b -C 1, then from Eq. 
effective resistance is thus 

ReS ~ 2’ R ^ b. (3.13) 

This equation does not assume 612* <C 1 nor Rfl* 1. 


(3.5), we have V = r2*T* = ^'i/>oQ^,{R/by. The 


4 Black Hole with a Disk 


In their original paper [3], Blandford and Znajek derived the stream equation in the Kerr spacetime 
[the analog of Eq. (2.4)] and laid out a perturbative method for solving it. One begins by expanding 


the flux function, current, and angular velocitjj^in the spin of the black hole as follows 

V. = V(r.») + o(A 


(4.1) 

(4.2) 

(4.3) 


where X, Y, and W are independent of the spin a|^ [We work in Boyer-Lindquist coordinates 
{t,r, 6, (/>).] Since the current and angular velocity vanish at a = 0, the leading-order flux function 


Kerr, or more generally any stationary axisymmetric (circular) spacetime, these quantities are invariantly 
defined in the same way as in flat spacetime m- 

^®The leading scalings follow from the fact that no energy can be extracted from a nonrotating black hole. [Note 
that Eq. ( 2.3[ ) holds in Kerr.] The scaling of the error terms may be seen from the Znajek condition (e.g., Eq. (104) 
of Ref. IT 
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Our notation differs slightly from BZ who use = tjj — I = X + 0{a/My 


13 









X must be a vacuum Maxwell solution in the Schwarzschild spacetime. At next order 0(a), the 
equations are automatically satished by the form of Eqs. (4.1)-(4.3)—the entire content of the 


force-free equations is the statement that Y and W are functions of X alone. The Kerr force-free 
stream equation appears hrst at 0{a^), which we do not consider in this paper. 

To hnd force-free solutions to 0{a), then, one needs a vacuum Schwarzschild solution X as well as 
two further equations to £x Y and W. The hrst comes from a universal relationship that always 
holds on the black hole horizon, the so-called Znajek condition 


Y = [W 


- ) sin 6* deX 


at r = rn, 


(4.4) 


which we present here to the relevant order in the spin. This condition guarantees that the helds 
are regular at the horizon and also follows from the stream equation (e.g.. Ref. m)- Note that | 
is the value of W corresponding to the horizon angular frequency. The second equation, according 
to a key idea of the BZ approach, is to be determined from an analogous force-free solution in flat 
spacetime. If the goal is simply to satisfy the force-free equations to 0(a), then any postulated 
second relationship between X,Y, and W will do. But this is an extremely large freedom and may 
correspond to spurious linearized solutions that arise from no family of exact solutions. In order 
to eliminate this freedom, BZ postulated that a second relationship should be determined from an 
exact hat spacetime solution that agrees with the linearized solution at large r, where the Kerr 
spacetime becomes approximately hat. In our case, this second relationship will be Eq. (2.18) with 
the lower sign appropriate to outgoing hux in the northern hemisphere. 


Y = -WsmOdoX 


as r —)■ oo. 


(4.6) 


Here and below, we work in the northern hemisphere 0 < 0 < 7r/2, with the southern helds 
determined by rehection. 


4.1 Flux Function X 


The hrst step of the BZ method is to hnd a vacuum (J = = 0) hux function X in Schwarzschild 
that corresponds to the “same” physical conhguration as the hat spacetime hyperbolic force-free 
solution ip evaluated at / = = 0. In the cases treated by BZ, there was an obvious “same” 

solution, but here we must exercise more care. In general, there is no canonical mapping between 
held configurations on diherent spacetimes, but in this case, we have just enough structure to 
make one. Using the facts that (i) the hux functions must satisfy the vacuum Maxwell equation in 
their respective spacetimes (hat and Schwarzschild) and (ii) these spacetimes both exhibit spherical 
symmetry, we can make a unique identihcation mode by mode in an expansion in angular harmonics. 
In Appendix we follow this general approach and derive a family of solutions parametrized by 
the Schwarzschild coordinate radius ho of the disk, expressed in series expansions valid for r < bo 


and r > bo- The hux function for this family is given by Eqs. (B.22), (B.24), and (B.25). 
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For the discussion here, we need only the following properties of the solution: 

(4.6) 

(4.7) 


X{r = th, 0) ^ Xh sin^ 6, 

X(r —)■ oo, 6*) = Xo(l — cos6*). 


These indicate that the held is approximately vertical on the black hole and approximately radial 


(monopolar) at inhnity. We have multiplied the solution (B.22) by an overall constant Xq repre¬ 


senting the magnetic hux in one hemisphere. Equation (4.7) follows from (B.22b) as r —)■ cx), noting 
that Rf{r) ~ r~^. Equation (4.6) follows from (B.22a), keeping only the leading k = 1 term, noting 
that Rf{r) = r^. Dropping the higher-order terms is an excellent approximation provided that the 
disk does not extend too close to the black hole. For example, for bo > 6M, the subleading term is 
at least a factor of 10“^ smaller. The constant Xh is the total hux on one hemisphere of the black 
hole horizon and is given by 


Xm = CX, 


8M2 


07 


bo-M+ ^bo{bo-2M) 


2 ’ 


where C{M,bo) is a constant of order unity given in Eq. (B.25) [see also Eqs. (B.39)-(B.40) 
bo = 6M (the ISCO radius), C ^ .735, while for bo ^ M (disk very far), C = 1. Thus, 


Xo 


bo=GM 


Xh 

.060, 

-^0 


2M^ 


(4.8) 


For 


(4.9) 


At bo = 6M, the bo ^ oo expression gives Xh/Xq ^ 0.056, quite close to the true value of 0.060. 


4.2 Current Y and Angular Velocity W 


Plugging Eq. (4.6) into (4.4) gives a condition at the horizon, while plugging Eq. (4.7) into (4.5) 
gives a condition at inhnity. These, respectively, are 


Y = 2Xh {W — siiY 9 cos9 atr = rH- (4-10) 

F = —XoVF sin^ 0 as r —)■ cx). (4.11) 


These equations hold at diherent radial coordinates and cannot yet be compared. However, we may 
use Eqs. (4.6) and (|4.7[) to eliminate 9 in favor of X, 


y = 2 (H7 - 1 ) .Y Vl - x/x„ (from the horizon; X < Xh), (4-12) 

Y = -WX{2 - X/Xo) (from inhnity; all X). (4.13) 

Since Y and W are functions only of X, the equations may now be combined, giving 


F(X) 

W{X) 


-WX{2-X/Xo), 

1 ^l-XjXn 

4l-X/(2Xo) + y/l-X/XH’ 


(4.14) 

(4.15) 
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which hold on the overlapping domain of validity, X < Xh- For X > X^, 'we take F = hF = 0, as 
required by reflection symmetry and smoothness. Note that this does not entail any current sheet 
since Y = W = aX X = Xh from Eqs. (4.14) and (4.15). Recall that 1/4 is the value of W 


corresponding to the horizon angular frequency. 


A similar procedure was performed previously by Ref. m (see also Ref. 

Eqs. (4.14) and (4.15) in the distant-disk limit bo ^ M (equivalently Xh Xo 
the terms with X/Xq disappear, since X < Xh in those expressions. They postulated Eq. (4.13), 


who arrived at 
In this limitm=. 


which we derived here from the corresponding flat force-free solution (2.8) according to the BZ 


method. They found the vacuum flat solution ( |B.18[ ) in generality and used it to hnd the near¬ 
horizon behavior of the Schwarzschild vacuum solution when ho S> M. (In this limit, (7 = 1 and 
b = bo, and the many subtleties discussed in Appendix]^ do not arise.) We have found the vacuum 
solution everywhere outside the horizon and for arbitrary bo , an effort which reveals that the bo ^ oo 


expression in Eq. (4.9) in fact works well even at the ISCO radius bo = 6M. Equations (4.14) and 


(4.15) without the XjX^ terms were also obtained in Ref. [3S] by considering strictly vertical held 
lines and imposing an integrability condition for the second-order perturbation. 


4.3 Black Hole Jet Properties 

As in the exact hat solution, the approximate Kerr solution produces a conical jet. The boundary 


of the jet is the last held line touching the black hole, labeled by X = Xh- From Eq. (4.7), this has 
an asymptotic 6*-angle of 


^j’et = arccos ( ^ ^ ~ 


bo = 6M, 

bo > M. 


(4.16) 


At large r, the power per unit 6 in the jet is 


dP x^n% 


where 


w{9) = 


y^cos 9 — cos 9*^^ 


y/1 — cos cos^ I -|- y^cos 9 — cos 6*j* ^ ’ 
and VLh = a/(4M^) at this order of approximation. 


(4.17) 


(4.18) 


The total power may also be computed from Eq. (2.3), which holds in any stationary, axisymmetric. 


circular spacetime m- We can compute numerically for all hnite values of bo or get an analytic 
expression when 6o 3> M (equivalently Xh Xq), 


7 X 10-5 X2 Qjj ^ 


.02Xjjnl 


P* 

Fet 


TT 


f^Xn 


X ( ^ dx = 

yi+^i-x/XH J 


17-6 log 16 ^2 o 2 


67r 




bo = QM, 

bo > M, 


(4.19) 


16 
























where we have used Eq. (4.9). Note that (17 — 61ogl6)/(67r) .02, so that the large-6o expression 

in fact gives an excellent approximation for the bo = QM answer. These expressions for the power 
include only one jet; the total power is a factor of 2 larger. 


4.4 Circuit Model 

The voltage drop across the jet is the integral of Qd'ip. We can compute this analytically in the 
limit 6o 3> M (equivalently Xh -C Xq), 

V = [ndij = nH [ 5 o > m . ( 4 . 20 ) 

J Jo l+^Jl—X/XH 

In this limit, the net outward current —I reaches a maximum value of \Xh^h at X = \Xh- We 
regard X < \Xh as the “bulk” flow and X > \Xh as the return flow and write 

= (4.21) 

This is the current in the circuit, and hence the effective resistance is 

= 2(log4 - 1) ^ .77. (4.22) 

Ibulk I 

This analysis holds for 6o S> M but provides a good approximation even for bo = 6M. 


5 Stars vs. Black Holes 


We conclude with a side-by-side comparison of the star (Sec. and black hole (Sec. solutions. 
The main qualitative difference is that the star’s jet is surrounded by a current sheet, while the 
black hole’s jet is not. To compare more quantitatively, we must specify how the parameters are to 
be chosen in each case, make the same approximations in both compact object (CO) solutions, and 
compare only invariantly dehned quantities. On the one hand, the star solution is parametrized 
by radius R, disk edge b, and angular velocity O*. On the other hand, the black hole solution is 
parametrized by black hole mass M (also proportional to its size), Schwarzschild coordinate disk 
edge bo I and angular velocity VLh- We will take Oco = = O* and assume that the disk edge is 

far from the compact object, b ^ R and bo ^ M; then, we may take bco = b = bo without any 
coordinate ambiguity. There is no unambiguous way to compare the stellar radius R with the black 
hole mass M. Instead, we will demand that the opening angles of the jets are the same. From 
Eqs. (3.6) and (4.16), we have = R/b and = 2M/bo. Hence, setting them equal gives 




R = 2M, 


(5.1) 
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Figure 4: Comparison of jets from stars and black holes immersed in the magnetic held of a thin 
disk. On the left, held line angular velocity as a function of magnetic hux, normalized by the 
compact object velocity Oco and the magnetic hux on the compact object tl/co, respectively. In 
the star case, the held angular velocity equals that of the star, while in the black hole case, it is 
variable, peaking at half the black hole frequency on the pole and vanishing on the equator. The 
star current and held velocity do not vanish at the edge of the jet, indicating the presence of a 
current sheet there. To make a direct comparison, we have assumed that the compact objects are 
slowly spinning and that the the disk is far away. More generally, the current in the star’s jet is a 
nonlinear function of magnetic hux, as indicated in Eq. (2.8). 


which happens to be the Schwarzschild coordinate radius of the event horizon. This also makes 
the total magnetic hux Tco = ^* = Xh = ^'ipo{R/by on the objects equal—see Eqs. (3.5) and 
(4.9)—while identifying Xq = 'i/jq. With these identihcations, we have 


P* = 

Met 


(Tcof^co)^ _ 17-6 log 16 2 

-, Met- - PZ - (^CO“COj 


2^ ^ -.e. ^ i^/et-lOxP;, (5.2) 

The total power of the star is approximately ten times larger than the black hole’s power. In this 
limit, we have the simple expressions (3.12) and (4.22) for the ehective resistance of the circuit, 

5, P:ff = 2(log4-l)^.77. (5.3) 


r:s = 


In analytical estimates, one sometimes assigns the impedance of free space Pfree = 1 (~ 377 Ohms) 
to the ambient space and to the black hole. This would suggest = 1 and = 2, which is oh 
by a factor of about 2. However, the basic intuition that the black hole provides extra resistance is 
supported by the exact solution. 

Since the black hole solution is valid only for MQh <M 1, we should similarly assume PO* <M 1 
in the star case. With all these approximations. Fig. plots the angular velocity and current 
functions of the star and the black hole. The physical setup is nearly identical, and the compact 
objects are ehectively vanishingly small in the approximations taken. Nevertheless, the distant jet 
has dramatically diherent properties in each case. In this way, the force-free magnetosphere carries 
information about the nature of the central object to observers located at inhnity. 
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A Flux Function for Constant Q 


We now evaluate the integral in Eq. (2.8) for i/j in the case that r2(n) = H is a nonzero constant. 
We must then require that |6r2| G (0,1) and hnd 




V’o 


£ 


^bVL {hVL + 1) \ 


arcsm 


bVLv? + 1 


’612 + 1 


+ D, 


(A.l) 


where £{(j),m) is the elliptic integral of the first kind. The integration constant D must be fixed 
such that V’ vanishes on the pole, "0(0) = 0- This gives 


'ip{u) = 




y/bQ {bQ + 1) 




arcsm 


bO.u'^ + 1 


612 + 1 




4 ’ 612 + 1 


(A.2) 


B Vacuum Solution in Schwarzschild 


In this seven-part appendix we discuss general procedures for identifying vacuum (source-free) 
Maxwell solutions in flat spacetime with corresponding solutions in Schwarzschild spacetime and 
present an explicit construction for the hyperbolic case of relevance to this paper. 


B.l Overview 


In the vacuum case (J = 12 = 0), the stream equation (2.4) in flat spacetime reduces to 

+ dl'ijj — cot 9 deip = 0, 


(B.l) 


where now we work in spherical coordinates (r, 6^). The corresponding equation in Schwarzschild 
spacetime (with Schwarzschild coordinates) is 


r(r — 2M) + 2M drip + dlip — cot 9 dei> = 0, 


(B.2) 


where M is the black hole massj^ These equations are separable, have the same angular part at all 
r, and agree entirely as r —>■ oo. These facts allow one to make a canonical identification between 
solutions in the M = 0 and M ^ 0 cases, as we now explain. 


^^Equations (B.l I and (B.21 are equivalent to the vacuum Maxwell equations, in their respective spacetimes, for 
an axisymmetric vector potential in a gauge where only = pj is nonvanishing. They differ from the equations 
for the magnetostatic scalar potential used in Ref. [39] and many other references. 
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A physical flux function must vanish at both poles. As discussed in Appendix B.2 below, a 


sufficiently regular solution of (B.2) vanishing at both poles can be expanded, at least locally in r. 


as 


^ [Ci'Rfir) + CfR^r)] 0,(9), 


(B.3) 


£=1 


where Cf 


and Cf 


are (real) constants. The angular harmonics 0£ vanish at the poles, while Rf is 
regular at the horizon (but not at infinity) and Rf is regular at infinity (but not at the horizon). 
Solutions of the M = 0 equation (B.l) admit an expansion analogous to that in Eq. (B.3), with 
identical angular harmonics Qe{0). The radial harmonics are different but may be identified with 
corresponding M f 0 radial harmonics by demanding that they agree as r —)■ oo, where the 
equations agree. Thus, one may associate solutions of Eqs. (B.l) and (B.2) by expressing them in 
the form (B.3) and demanding that the coefficients Cf be the samej^ 

While general, this procedure is not particularly convenient for force-free magnetospheres. Physical 
magnetosphere models typically contain a current sheet, which appears as a cusp in the flux function. 
(For example, the split monopole has -0 = 1 — |cos 6^1.) While such functions can be expanded in the 
series (B.3), the representation is not intuitive, and the cusp limits the convergence. Furthermore, 
computing the coefficients Cf involves evaluating integrals of the form (B.9), a difficult task in 
general. To overcome these difficulties, we will use a method appropriate for magnetospheres with 
equatorial reflection symmetry, where we work only on the northern hemisphere 9 G [0,7r/2], with 
the southern fields to be determined by reflection. This frees us to consider functions that vanish 
only on one pole, which means the £ = 0 harmonic 0o(^) = 1 — cos 6* is now admissible. The 
corresponding radial harmonic is Ro{r) = 1, with the other linearly independent solution divergent 
at both the horizon and at infinity. 


The method takes advantage of the fact that the flat space radial harmonics are pure powers of 
r, Rf = and Rf = r~^, making Taylor expansion equivalent to mode expansion. We Taylor 
expand the flat solution about r = 0 and r = cx) to produce expansions valid for r < b and 
r > b, respectively. We then promote each expansion to a solution in Schwarzschild by the above 
procedure. The series converge on either side of a Schwarzschild coordinate radius r = bo{M,b) 
that we compute explicitly. The overall normalization of the outer solution is fixed by demanding 
that the flat and Schwarzschild solutions agree asymptotically where both spacetimes are flat, while 
the overall normalization of the inner expansion is fixed by demanding continuity at r = 6o- The 
remainder of this appendix follows this procedure. 


B.2 Angular and Radial Harmonics 


The separation of Eq. (B.2) has been discussed previously by Refs. |2S1 SDl SI]- Our treatment 
establishes notation for this paper, provides some additional details, and finds the normalization of 


^®While this can always be implemented in a sufficiently small (but finite) domain, in typical problems, Cf and 
Cf will be chosen differently on different domains (e.g., matched at some radius). In this case, there remains the 
ambiguity of how to choose the domains. This ambiguity does not arise in the related method we use in this appendix. 
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the radial harmonics required to properly reduce to the canonical normalization in flat spacetime. 
The product solutions Ri{r)Qi{9) of Eq. (B.2) satisfy 


d f 1 dOA iji + l) ^ 

dO \sin0 dO ) sin 6* 
d 
dr 


2,M \ dRi 
dr 




(B.4) 

(B.5) 


where we name the separation constant i{i + 1). Since the angular equation (B.4) is of Sturm- 
Liouville form with weight esc 6, we can expeetp^ a discrete, infinite set of eigenvalues whose eigen¬ 
functions are complete and orthogonal with respect to esc 6 dO. We And that the solutions vanishing 
at both poles are given (in a convenient normalization) by the hypergeometric function, 


02fc-l(6') 

02fc(^) 


2-^1 

2-^1 


—k, k — 
—k, k -\- 




1 


COS 6 , 


(B.6) 

(B.7) 


where k (and hence t) is a positive integer. The 0^ are proportional to the Gegenbauer polynomials 


0,(0) = 


2i"l 


2i"l 


+ 11 2 n 

-; -;cos 6 


2 2 ' 


+ 13 2 a 

; -; cos 6 




2v/7r 


2 ’ 2 ’ 2 


COS0 = 




£ odd, 

(B.8) 
^ even. 


From the known orthogonality properties of those polynomials, one can directly check that the 
0, are orthogonal with the expected weight csc0d0. Given a sufficiently regular function T(0) 
vanishing at the poles, a series representation can now be computed as 


£=1 


‘ /„'e((9)^csc9d9 


(B.9) 


We will also And it useful to consider the £ = 0 solution, which is a -1- /3 cos 6 for constants a and 
(3. We can make this function vanish at one pole and hence use it in a single hemisphere. We work 
with the choice 


0o(0) = 1 — cos 9. 


(B.IO) 


This function is not orthogonal to the functions 0,>i defined in Eqs. (B.6) and (B.7), and we do 
not regard it as part of that set, which is expected to be complete only for functions vanishing at 
both poles. 
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The function esc 0 = 1/ sin 6 is not sufficiently regular for the standard theorems to directly apply. 
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In a convenient normalization, the solntions of the radial eqnation (B.5) regnlar at the horizon and 
inhnity, respectively, are given for ^ > 1 by 


Rfir) = -(-2M) 


RUr) = 


T{2i+1) 

r(£+|) 


>Fi 


+ 2,l-£;3; 


2 M 




iFi 


+ 2 ,£; 1 ; 1 - 


2 M 


log 1 - 


2 M 


+ Pi 


2M 


where the are the polynomials recnrsively dehned by 


X 


Pi{x) = a; + -, 


Poix) = Ax'^ — x^ 


X 


6 ’ 


PAx) = 


{2i - - l)(2a; - 1) - l]a:P^_i(a:) 2)x‘^P^_2{x) 


(£+!)(£- 1)2 


Note that the leading term in Pi is 


Fix) — ^ 2) 2e I nl^22-i\ 


(B.ll) 

(B.12) 

(B,13) 

(B.14) 

(B,16) 

(B.16) 


The normalization of Rf is chosen so that as r —)■ cxd hxing M (or equivalently M —)■ 0 hxing r), 
we reduce to the simple form 


Rf{r) = r 


.i+i 


R> (r) = r" 


as 


r —)■ 00 or M —)• 0 . 


(B.17) 


This makes the M 7 ^ 0 normalization of Eq. (B.ll) the “same” (in the sense discussed in Sec. B.l) 
as the natural M = 0 normalization Rf{r) = and Rf{r) = r~^. When M 7 ^ 0, Rf{r) diverges 
logarithmically on the event horizon, r = 2M. 

Note that the M = 0 eigenfunctions in Eq. (B.17) differ from the usual A and that arise in 


solutions of the Laplace equation. The reason is that we work with the flux function rather than 
the magnetostatic scalar potential (see discussion in footnote mi- 

For £ = 0, the general solution is a + /3[r + 2Mlog(r — 2M)] for some constants a and £3. Unlike 
the f' > 1 cases, the £ = 0 solution is either everywhere regular (/3 = 0 ) or divergent at both the 
horizon and inhnity (/d 7 ^ 0 ). 


B.3 Mode Expansion of the Flat Solution 


Setting '^0 = 1 for simplicity, the hat vacuum solution ip, Eq. (2.9), is given in spherical coordinates 
by 


A{r, 9) = 1 


\ 


r \ 2 

b) 


+ ''i 


+(()hos^». (B,18) 
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In principle, one could proceed by integrating tjj against the angular harmonics 0^ from Eq. (B.9) 


to determine the coefficients Cf in Eq. (B.3), but this approach is intractable. We will take an 
alternative approach based on the observation that the radial eigenfunctions in flat spacetime are 
and r~^ exactly, and therefore correspond to terms in a Taylor expansion around r = 0 and 
r = oo, respectively. Thus, by Taylor expanding, we can expect each term to come with an angular 


function 0^ and a coefficient that can easily be read off. For example, expanding (B.18) near r = 0, 
we hnd 

sin^6* /r\2 sin^ 6* , /r\4 sin^6* , ^ /r\6 

-(3 + 5cos20) (15 + 28 cos 26* + 21 cos 40) +..., 

(B.19) 




each term of which has 0-dependence proportional to an angular function Qi. The full series is 


= 


E 


r(k-i) ,r 


2k 


02fc-l(^), 


r < b. 


^ r(A; + 1) \b. 

We may similarly expand around r = oo and read off the complete Taylor series as 


'ip{r, 0) = 1 — cos 0 — ^ 


2k 


T(2k+l)(h 


r > b. 


(B.20) 


(B.21) 


“ r(fc -i-1)2 

Here and hereafter, we work in the northern hemisphere 0 < 0 < 7r/2. Notice that the near series 


(B.20) involves only odd multipoles i = 2k + 1, while the far series (B.21) involves even multipoles 


i = 2k. (Note that 1 — cos0 is the 6 = 0 harmonic.) These series do not match term by term and 
must be summed for comparison at r = 6. 

We have established the radii of convergence of these series by using the ratio test at the special 
value of 0 = 7r/2 and conhrming numerically that the radius is 0-independent. This procedure is 


the exact flat solution (B.18). 


described in Sec. B.5 below. Equations (B.20) and (B.21) provide a dual series representation of 


B.4 Schwarzschild Solution 

To promote the flat solution to Schwarzschild spacetime, we send —>■ Rf{r) and r~^ — Rf{r) 


in Eqs. (B.20) and (B.21) and also allow for a relative normalization C, 

2k 




0 ) = < 


r(fc + 1 ) \b 

1 — cos 


^r(fc + i)2 


In Sec. B.6 below, we determine the radius of convergence bo to be 


bo = b + M+—. 

Ab 


r < bo, 


« - i; ((!) r > 


(B.22a) 

(B.22b) 

(B.23) 
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This indicates that the disk terminates at the Schwarzschild coordinate r = bo- That bo ^ b 
should not be surprising, since there is no canonical identihcation between coordinates on different 
manifolds. (Note, however, that bo ^ b as M ^ 0, as required.) Inverting Eq. (B.23) gives 


b = 


bo — M + \/bo{bo — 2iVf) 


(B.24) 


We regard bo as a free parameter in the Schwarzschild solution. 

In principle, the normalization C may be determined by summing both series at r = 6o- This is 
tractable at the special value oi 9 = ti/2 where the angular harmonics simplify. The even harmonics 
vanish, showing that >b) = 1. The odd harmonics are equal to unity, so we have 


1 

C 




r(i-i) 


2a/T T{k + 1) 


2k 




(B.25) 


This series converges excruciatingly slowly, making direct computation with Eq. (B.25) impractical. 
Its terms only fall off like so it is necessary to sum the first ~ 2, 600 terms to get to within 1% 


accuracy. In Sec. B.7 we derive an equivalent, more practical representation using series acceleration 
techniques, Eq. (B.39). 


To complete the proof that Eq. (B.22) is a smooth solution in Schwarzschild one should in principle 


check that the series match at all angles 6 using the expression (B.25) for C determined at 0 = 7r/2 


We have verified this numerically in broad brush, but it is difficult to be conclusive because of the 
slow convergence of at least one of the series at every value of 6. One can expect the series to match 
because of the unique, natural association between solutions in Schwarzschild and solutions in flat 
spacetime discussed in Sec. |B.1[ This association acts over finite neighborhoods and up to overall 


normalization. Thus, we can expect only a relative normalization when we promote our single flat 
solution in the two different neighborhoods. 


B.5 Flat Radius of Convergence 


We now discuss the radius of convergence of the series representations (B.20) and (B.21) of the flat 
solution. We will need the asymptotic approximation for the hypergeometric function given on p. 
77 of Ref. 


2F1 


a + A, 6 — A; c; 


1-z 

2 


A-s-oo r(i & + A)r(c) 
r(i)r(c-6 + A) 

X A-1/2 


{l+e-T 

j,-(A+.)«] [l + c>(|A-‘|)], 


a-b-1/2 


(B.26) 


where the upper or lower sign is taken according to O' 2 ; ^ 0, while ^ is defined hy = z-\- ^z^ — 1. 
(Note that = 2 ; — ^z^ — 1.) For real arguments 2 ; = a + Oi, one considers 2 ; = a ± ie and lets 
e —)■ 0, where the result is independent of the choice of ±. 
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This formula allows us to find the asymptotic behavior as fc —>■ oo of (B.6) by setting a = —1/2, 
b = 0, c = 1/2, X = k, and z = 1 — 2 cos^ 9, resulting in 


02fc-l(6') 


fc^oo T{k + 1) 1- e 


o2ie r 


r(fc+l)V 8A: 

Hence, we may rewrite the sum ( |B.20[ ) in the asymptotic form 

fc—>oo rf/c 2) /r 






(B.27) 




2k 


1 — e' 
327r/c 


2i6» r 




-2i6»\-^+l/2 


(B.28) 


This expression is real at all 9 but involves a delicate cancelation of complex phases at generic values 
of 9. The ratio a^+i/afe in general oscillates wildly on account of these phases, and the ratio test 
only gives weak bounds on the domain of convergence. However, at the special value oi 9 = tt/2, 


L = lim 

Ofc+l 

= lim 

k^oo 

Oik 

k—^oo 


2k 


k 


T’ \ 2 

2k + lS k + l\h 


r\ 2 

b. 


(B.29) 


Thus, the ratio test shows that the series converges for all r < b. The test is inconclusive for r = b, 
but the technical convergence of the series there is unimportant for our purposes, since the value 
can always be computed by working at r = 6 — e and letting e —?• 0. For this reason, we write r < b 


as the domain of validity of the series (B.20). We have not proven that the radius is r = 6 for all 9, 


but this is straightforward to verify by numerical experiment. Similarly, one may check numerically 
that the outer series converges for r > b. 


B.6 Schwarzschild Radius of Convergence 


Finding the radii of convergence of the Schwarzschild series (B.22) requires the asymptotics of the 


radial functions (|B.11|), which also follows from Eq. (|B.26|). For the inner series, we take a = 1, 

r — •\/r(r — 2M) 


6 = 2, c = 3, A = 2k, and z = 1 — r/M to find 
/c—>00 /3^F(2fc-l) 


^2k-l^ 


Tik F(2fc + 1) (1 — \/l — e 

Hence, we may rewrite the sum as 


3 ^ = 1 


M 


(B.30) 


i^{r,9) = Cj2h, 


k=l 


2a/F T{k + 1) \b 


2k 




fl2Vi(r)e2t-i(9), 


(B.31) 


and use Eqs. (|B.27|) and (|B.30|) to obtain the asymptotic behavior of bk, 

k^{2k- 1)2 (1 - e€) 


7 k^oo 
Oh ~ 


27r V6/ 


T{2k + 1)2 /2M\ / 1 - e2*« (e 

F(4fc-1) V~r; V1 - e-2€ A;2(2fc- 1)2 


-2kiB ^g(2fc-l)i6»j 
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(B.32) 
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As in the asymptotic formula (B.27) for the angular functions, the ratio test gives only weak bounds 


for generic values of 6, on account of complex phases. At the special value of 0 = 7r/2, however, we 
have 


L = lim 

k^oo 


^fc+l 


/My 

r — ■\/r{r — 2M) 

bk 

[2b) 1 1 


^ M 


(B.33) 


where the last step follows from Eq. (B.30). Demanding L = 1 gives the formula (B.23) for bo, 


where L < 1 for all r < bo so that convergence is assured there. In the same spirit as the discussion 


below Eq. (B.29), we regard r < bo as the domain of validity. Similarly, we expect r >bo to be the 


domain of validity of the outer series, which is straightforward to conhrm numerically. 

B.7 Series Acceleration for Matching Coefficient 


To facilitate numerical evaluation of the normalization C in Eq. (B.22a), we partially resum the 

rgent expression. 

1 T{k-^) /!' 


series (B.25) to produce a more rapidly convergent expression. The series may be written as 

2k 


1 

c 




— 


k=l 


2^/n r{k + 1) \b 




(B.34) 


asymptotic behavior (B.32), we have 


The Ck are equal to the b^ dehned in Eq. (B.31) evaluated at r = 6o and 6 = 7r/2. From the 

2^^T{k-\)T{k+\)T{2k-l) 


k—>oo 1 

Ck ~ 




71y/2b-M {2b + 
where we used the fact that at r = 6o, 

bo - \/bo{bo - 2M) 


y/kT{4:k-l) 


(B.35) 


= 1 


M 


M 


(B.36) 


with b a function of bo as in Eq. (|B.24|). This large-fc asymptotic behavior is not yet exactly 

summable, so 
simplify it to 


summable, so we appeal to Stirling’s approximation r(^) in order to further 


k^OO 1 

Ck ~ dk, 


du — 


hhl 


V2b-M {2b + ' 


(B.37) 


We now recognize the weak k falloff as the reason behind the slow convergence of the series. 
But the sum of k~^ is the Riemann zeta function C(^)) so we have 


k=l 


^ _4 _ b^ _ / 3 

^ ~ V^^j2b-M {2b + V2 


(B.38) 
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Thus, we may instead write 
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C 


4 bhl /3\ 

A V26 - M (2b + 


oo 

+~ Cfc) . 
A:=l 


(B.39) 


The series of correction terms converges quite quickly. If we keep only the hrst term k = 1, the 
estimate of C is 


1 

C 


+J_ 

2b^ \/^V2b-M (2b + 



(B.40) 


This truncated expression is generally accurate to within about 5%. For example, if the disk extends 
to the ISCO, bo = risco = 6M, we get an exact value of C~^ = 1.3603... from Eq. (B.39) and an 
approximate value of C~^ ~ 1.289 from Eq. (B.40). For b ^ M, we know that C = 1 exactly, and 
Eq. (|B.40l) gives C~^ .955. S 
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